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Interfaces in oxide heterostructures always provide a fertile ground for emergent properties. 
Charge transfer from a high energy band to a low energy opponent is naturally expected, as oc¬ 
curring in semiconductor p-n junctions. In this study, several exceptional physical phenomena have 
been predicted in (YFe 03 )n/(YTi 03 )TO snperlattices. First, the charge transfer between these Mott 
insulators is in opposite to the intuitive band alignment scenario. Second, hybrid ferroelectricity 
with a moderate polarization is generated in the n = 2 magnetic superlattice. Furthermore, the 
ferroelectric-type distortion can persist even if the {ABOs) 2 /{AB' 0^)2 system turns to be metallic, 
rending possible metallic ferroelectricity. 

PACS numbers: 73.20.At, 75.70.Cn, 77.55.Nv 

I. INTRODUCTION II. MODEL & METHOD 


Oxide heterostructures provide a unique and fertile 
ground to study emergent physics of correlated electrons 
and a promising route to design new devices using quan¬ 
tum effects. At the interfaces between oxides, collec¬ 
tive behaviors of electrons can be contrastively different 
from their original roles in parent materialsii^. For ex¬ 
ample, previous studies have revealed plethoric phenom¬ 
ena relevant to electronic reconstructions, e.g. metal- 
insulator transition and enhanced Neel temperatures 
in LaMnOa/SrMnOa snperlattices^^— , two-dimensional 
electronic gas and superconductivity in LaAlOa/SrTiOs 
heterostructures2i^, orientation-dependent magnetism in 
LaNiOa/LaMnOgSid^ and LaFeOa/LaCrO j^^d^ super lat¬ 
tices. 

Charge transfer is one important driving force for elec¬ 
tronic reconstructions, which tunes the local electron 
density as well as the interfacial electrostatic field. In 
the light of the band alignment scenario, the heights of 
Fermi levels are the decisive factor for charge transfer. 
Naturally, to reach a uniform chemical potential across 
the interface, electrons will leak from the high energy 
band to the low energy opponent, which is well known in 
semiconductor p-n junctions. 

Different from the nearly-free electrons, the bands 
of correlated electrons are not strictly rigid but some¬ 
what “soft”— which may lead to emergent phenom¬ 
ena beyond the simple band theory. In this work, the 
(YFe03)„/(YTi03)„ snperlattices have been studied us¬ 
ing the density functional theory (DFT). Several un¬ 
expected physical phenomena, e.g. the charge transfer 
against the band alignment scenario and hybrid mul- 
tiferroicity, have been predicted. In addition, due to 
the origin of improper component of polarization, the 
ferroelectric-type distortion in {AB 03 ) 2 l{AB' 0^)2 su¬ 
perlattices can persist even if the system turns to be 
metallic. 


Both YFe03 and YTiOs are Mott insulator a^^d^ . By 
neglecting the weak ferromagnetism due to tiny spin 
canting, the magnetic ground state of YFeOs is antifer¬ 
romagnetic (AFM) ^’^d^ . In contrast, YTi03 is ferromag¬ 
netic (FM)i^. These two materials share the common 
A-site cation Y^+, as well as the identical space group 
of crystal structure (orthorhombic Pbnm, see Fig.[T](a)). 
Their lattice constants (a, b, c) in unit A are: (5.283, 
5.592, 7.603) for YFeOai^ and (5.338, 5.690, 7.613) for 
YTiOg^. The proximate structures allow a high possi¬ 
bility for epitaxial growth of multilayers. In the follow¬ 
ing, the (YFe03)Ti/(YTi03)„ snperlattices are assumed 
to be grown on the mostly-used SrTi03 (001) substrate. 
To match the substrate, the in-plane lattice constants of 
YFe03 and YTi03 are hxed as 3.905 x ^/2 = 5.5225 A. 

All the following calculations are performed using the 
Vienna ab initio Simulation Package (VASP) ^^:^^ based 
on the generalized gradient approximation (GGA). The 
Hubbard UeS (= U — J) is imposed on Fe’s and Ti’s 
d orbitals using the Dudarev implementation^^. In the 
GGA-|-t7 calculation, the plane-wave cutoff is 550 eV. A 
7x7x5 Monkhorst-Pack fc-point mesh centered at F point 
is adopted for (YFe03)i/(YTi03)i and the parent mate¬ 
rials, while it is 6 X 6 X 2 for (YFe03)2/(YTi03)2. Both 
the lattice constant along the c-axis and inner atomic po¬ 
sitions are fully relaxed till the Hellman-Feynman forces 
are all below O.OI eV/A. 

For comparison, the hybrid functional calcula¬ 
tions based on the Heyd-Scuseria-Ernzerhof (HSE) 
exchange^ir— are also performed. Due to its extreme 
demand of CPU-time, the plane-wave cutoff is reduced 
to 400 eV. And the fc-point mesh is reduced to 3 x 3 x 2 
for (YPe03)i/(YTi03)i and the parent materials, while 
it is 3 X 3 X I for the n = 2 superlattice. 
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III. RESULTS & DISCUSSION 


A. Band alignment 

Before the study on superlattices, the parent materials 
are checked. According to the previous literatur o^^'^^ , 
proper C/es values C/pe = 4 eV and Uti = 3.2 eV are 
adopted in the following, if not noted explicitly. In our 
calculation, the magnetic ground states of YFeOa and 
YTiOs are G-type AFM (G-AFM) and FM, respectively. 
Both Fe^+ and are in the high-spin states. The 

relaxed lattice constants also agree with the experimental 
values quite wel l^^i^° . All these results guarantee the 
reliability of following calculations on superlattices. 

Then the substrate strain is imposed. The G-AFM 
order persists in strained YFeOa on SrTiOa substrate. 
However, for YTiOa, the strain can drive a magnetic 
transition to A-type antiferromagnet^^. The atomic pro¬ 
jected density of states (PDOS) of strained YFeOa and 
YTiOa are displayed in Fig. [TJb) andlljc), respectively. 
Both materials retain insulating with energy gaps of 
~ 2.4 eV for strained YFeOa and ~ 1.5 eV for strained 
YTiOa. 

According to Fig. ffib-c), the topmost valence band 
of YTiOa is from Ti’s one t 2 g orbital (whose position is 
denoted as /ixi), and the bottommost conducting band 
of YFeOa is formed by the spin-down orbitals of 
Fe (whose position is denoted as ^Pe). Both these two 
bands are very narrow, implying localized states. By 
aligning the deep energy bands of Y’s 4p and O’s 2s 
orbitals which should be identical in these two materi¬ 
als, a band alignment can be obtained straightforwardly. 
Ti’s occupied t 2 g band just locates within the forbidden 
gap of YFeOg, lower than the unoccupied conducting 
band of YFeOg. Mathematically, it can be expressed 
as /rxi < MFe- Thus an intuitive conclusion is that the 
charge transfer should be forbidden between these two 
materials, keeping the original Fe^“''-Ti^+ configuration 
across the interface. This band alignment is further con¬ 
firmed in the HSE calculations (Fig. [TKd-e)). 

Even if the above C/pe = 4 eV and t/xi = 3.2 eV are 
believed to be the most appropriate, it is interesting to go 
beyond the real materials by scanning the Hubbard Ues 
in a wider parameter space. By varying C/eff’s, the aligned 
band positions are summarized in Fig. [T](f). In general, 
/rxi decreases with increasing t/ex, while /rpe increases 
with increasing t/eff. This tendency can be understood 
based on the Hubbard model, since the occupied Ti’s t 2 g 
band is the lower-Hubbard band while the unoccupied 
Fe’s band is the upper-Hubbard one. Then, the intensity 
of charge transfer can be determined by comparing /ixi 
and /xpe, as summarized in Fig.[T](g). When both Gxi and 
C/pe are large, there is no charge transfer, as illustrated in 
Fig. HKb-c). In contrast, while in the small C/xi and C/pe 
limit, a partial or complete charge transfer should occur. 


B. Charge transfer 

Above analyses on charge transfer were based on 
the band alignment scenario. To verify this scenario, 
(YFe 03 )n/(YTi 03 )„ superlattices (n = I and 2) are 
studied. After relaxing the crystal structures, the mag¬ 
netic ground states of superlattices are determined. In 
both superlattices, the magnetic moment of Ti is (al¬ 
most) quenched, while the YFeOs layers retain the G- 
AFM configuration with a suppressed moment ~ 3.6 
/xb/F e. This magnetic reconstruction is due to the charge 
transfer, since Ti'^’*' is non-magnetic and the moment of 
high-spin Fe^’*' is 4 /xe/Fe. This scenario is further con¬ 
firmed by the PDOS’s (Fig. [Ha-b)). One of Fe’s upper 
Hubbard bands is occupied by one electron, while Ti’s 3d 
bands are fully empty, implying a complete charge trans¬ 
fer. Such a complete charge transfer is further confirmed 
using the HSE calculation (Fig. [2](c) and (d)). 

For physical comparison, the GGA-I-C/ calculations 
with lower C/eff’s have also been done. The charge trans¬ 
fer always occurs in superlattices’ calculation despite the 
value of C/eff’s. In other words, the prediction of charge 
transfer is not C/eff-sensitive. In addition, our GGA-I-C/ 
and HSE calculations give consistent results regarding 
the charge transfer, which are also in agreement with the 
previous LDA/LDA-I-C/ calculations and X-ray photoe¬ 
mission experiment on similar system LaTi 03 /LaFe 03 
superlattices^. Thus the prediction is not an artefact of 
the level of treatment of electronic correlations. 

Then how to understand such an unexpected charge 
transfer, which violates the band alignment scenario? 
Previously, similar charge transfer in LaTi 03 /LaFe 03 su¬ 
perlattices was attributed to the transition of Fe^"*" from 
the high-spin state to low-spin state^^. However, in our 
case, Fe^"*" remains in the high-spin state, giving ~ 4 /xb 
per Fe. In this sense, such a theoretical argument based 
on the low-spin state can not interpret the “anomalous” 
charge transfer predicted here. 

One may suspect that such an unexpected charge 
transfer is due to the reduced dimension of YFe 03 and 
YTiOs layers in superlattices. The band alignment 
shown in Fig. [IJb-c) is obtained according to the three- 
dimensional bulks, which may be reshaped into ultra-thin 
layers. Taking a tight-binding model for illustration, the 
reduced dimension can only tune the bandwidths but not 
the (central) positions of bands. According to Fig. [T](f), 
the shrinking of bandwidths can not reverse the band 
alignment when C/xi = 3.2 eV and U-pe = 4 eV. Thus, 
pure dimension reduction is not enough. 

Dimension reduction can also tune the Mottness by 
tuning the ratio of bandwidth and Hubbard U. Then the 
splitting between the upper and lower Hubbard bands 
can be shifted, which may alter the band alignment. 
However, this possibility can also be simply ruled out. 
The dimension reduction can only suppress the kinetic 
energy and thus prefers the Mottness, equivalent to in¬ 
crease Hubbard U. According to Fig. [Hf), the charge 
transfer is suppressed to zero with enhanced U. In short. 
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FIG. 1. (Color online) Properties of YTiOa and YFeOa. (a) Sketch of the common structure, (b-e) Projected density of states 
(PDOS) for YTiOs [(b) and (d)] and YFeOs [(c) and (e)]. The Fermi level in (b) is set as zero and the deep energy bands of 
Y’s 4p (around —23 to —22 eV) and O’s 2s orbitals (around —19 to —17.5 eV) are aligned for all four PDOS’s, as denoted by 
(green) dot lines. The Fermi level for each case is also marked by a (blue) broken line, (b-c) Obtained using GGA+17. (d-e) 
Obtained using HSE. (f) After the treatment of band alignment, the energy positions of near-Fermi-level bands of YTiOa (green 
upper-triangles, /rxi) and YFeOa (pink lower-triangles, /rpe) as functions of t/xi (upper horizontal axis) or l/pe (lower horizontal 
axis). In addition, the band edges of 1) topmost valence band of YFeOa (pink rhombs) and 2) bottommost conducting band of 
YTiOs (green circles) are also presented. Occupied/empty bands are marked by solid/open symbols, respectively, (g) Contour 
of charge transfer as a function of f/pe and t/xi according to the band alignment. The numerical label is on the top of box. 



Energy(eV) 


FIG. 2. (Color online) PDOS of superlattices, (a-b) GGA-l-17. 
(c-d) HSE. The Fermi level is positioned at zero (blue broken 
line). For n — 1, the deep Y’s 4p and O’s 2s bands are similar 
to their bulk’s correspondences, while for n = 2, clear split¬ 
tings are observed due to the electrostatic potential created 
by the charge transfer. 


the dimension reduction can not explain the unexpected 
charge transfer. 

All above evidences suggest that the reconstruction of 
Hubbard bands is the intrinsic origin for such an “anoma¬ 


lous” charge transfer. As illustrated in Fig. [2 for both 
Ti and Fe, the energy splittings between upper and lower 
Hubbard bands become much smaller than those in par¬ 
ent materials. The Fe^+, with half-filling 3d orbitals, is 
the most prominent candidate for the Hubbard repulsion 
between the spin-up and spin-down bands. However, the 
Hubbard repulsion will be suppressed when the electron 
density deviates from the half-filling, as in Fe^’*'. Similar 
mechanism works for Ti ions. In other words, the bands 
of YFeOa and YTiOa are not rigid neither but rather 
fragile against (virtual) hopping of electrons. Namely, 
although /rpe > MTi, both d^Fe/dn and d^Ti/dn are neg¬ 
ative, where dn denotes charge transfer from Ti to Fe. 

Similar mechanism should also work in 
LaFeOa/LaTiOa, which induces the charge transfer 
and low-spin state of Fe^+— . The low-spin state of Fe^’*' 
does not occur in our calculation due to the structural 
difference between LaFeOa and YFeOs. Following 
Ref. [ 2 ^ we can reproduce the transition from the 
high-spin state to low-spin state in LaFeOs by adding 
one more electron to Fe. However, the same treatment 
shows that the high-spin state of YFeOs is more stable. 
The highly distorted structure of YFeOs reduces the 
bandwidth of 3d orbitals, which prefers the Mottness 
and Hubbard splitting between spin-up and spin-down. 
In fact, the high-pressure experiment also found the 
high-spin state of orthoferrites is more stable when the 
A-site ion is smallii. 

Finally, it should be noted that although the Fe^'''-Ti'*^+ 
configuration also occurs in other systems, e.g. FeTiOs 
and Ti doped Fe 203 ^“— , the involved mechanism may 
be not simply identical. Although the electronegativity of 
different ions usually plays a crucial role^^, other factors. 
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especially the coordination environment and correlation, 
will also tune the charge transfer especially when the 
alignment between involved bands is subtle, as demon¬ 
strated in the present work. 


C. Hybrid ferroelectricity 

It was predicted that improper ferroelectric¬ 
ity could emerge in perovskite superlattices like 
(PbTi 03 )i/(SrTi 03 )A (LaFe 03 )„/(YFe 03 )r^, as 
well as the layer perovskites due to the 

modulation of non-polar antiferrodistortive modes. For 
these {ABOz)n/{-A-'BOz)^ the ferroelectric polarization 
appears to be nonzero when n is an odd number, while 
it is fully compensated between layers for even n’s^l. 

Our super lattices, in the type of {ABOz)n/{AB'Oz)n, 
is somewhat different. First, the crystal structures are 
analyzed. The space group of n = 1 superlattice is P^^jc 
, but that of n = 2 is Pmc2i. The point group of P^ijc 
is 2/m, and mm2 for Pmc2i. 2/m belongs to the non¬ 
polar point group while mm2 is a polar point group with 
its polar axis along the &-axis. Therefore, the n = 1 
superlattice should be non-ferroelectric while a finite po¬ 
larization is expectable in the n = 2 superlattice. Such a 
prediction based on symmetry analysis should be phys¬ 
ical robust and qualitatively independent on details of 
numerical calculations (e.g. 17eff)- 

The standard Berry phase method implemented in 
VASP is employed to evaluate the ferroelectric polariza¬ 
tion. As expected, the calculated polarization is zero 
for the n = 1 superlattice, but in the case of n = 2, a 
net polarization up to ^ 1.01 /iC/cm^ is obtained along 
the 6-axis. These results can be verified using piezoelec¬ 
tric force microscopy or optical second harmonic gener¬ 
ation, as done in the LaFe 03 /YFe 03 superlattices^^. In 
addition, the intuitive point charge model can also be 
employed to estimate the approximate value of polar¬ 
ization, which gives ^ 2.1 /rC/cm^, in agreement with 
the corresponding Berry phase values qualitatively. The 
quantitative difference is reasonable, since the Born ef¬ 
fective charge of ions can be different from their nominal 
valences. 

The ferroelectricity, together with the magnetic order¬ 
ing of Fe, makes the n = 2 YFe 03 /YTi 03 superlattice 
multiferroic. It is noteworthy that the previously stud¬ 
ied n = 2 LaFe 03 /LaTi 03 superlattice is non-magnetic 
due to the low-spin state of Fe^S. In this sense, the 
YFe 03 /YTi 03 superlattices can provide more physical 
functions. 

The origin of polarization can be visualized in Fig. |3l 
Due to the antiferrodistortive mode, both Y^+ and 
will move away from their corresponding high-symmetric 
positions in the high-temperature cubic structure. In the 
n = 1 superlattice, all displacements have their asymmet¬ 
ric opponents nearby (Fig. |3] (a)), thus can not generate 
a net polarization, as in the parent materials. In con¬ 
trast, in the n = 2 superlattice, the sequence of B-site 



FIG. 3. (Color online) Sketch of ferroelectric distortions. The 
arrows denote the displacements of Y®"*". {a) n — 1. The dis¬ 
placements are compensated between layers. The (b) positive 
and (c) negative ferroelectric distortion for n = 2. 


cations forms ...-Ti-Ti-Fe-Fe-... along the c-axis. This 
sequence is analogous to the ...-t-t-i-i"--- spin struc¬ 
ture in some type-II multiferroics (e.g. Ca 3 CoMnOg^, o- 
HoMnO s^^d^ , and BaFe 2 Se 3 ^), which owns parity. This 
parity, together with the antiferrodistortive mode, breaks 
the space inversion symmetry^. The displacements of 
Y^+ (and 0^“) sandwiched between the Ti-Ti (or Fe-Fe) 
bilayers and Ti-Fe bilayers are no longer asymmetric, giv¬ 
ing rise to a net in-plane polarization. Such distortions 
can be reversed by reversing the antiferrodistortive pat¬ 
tern, as illustrated in Fig. |3Kb) and[3Kc). Then the ferro¬ 
electric polarization is switched from positive to negative, 
as confirmed in the Berry phase calculation. 

Considering the similarity and difference, our result 
on {AB 03 )n/{AB^ 03 )n is complementary to the previ¬ 
ous studied (Ai?03)„/(A’i303)„, completing the theory 
of improper ferroelectricity in perovskite superlattices 
driven by the modulation of non-polar antiferrodistortive 
mode^l. In fact, the parity-related origin of ferroelec¬ 
tricity in {AB 03 )n/{AB' 03 )n superlattice is even more 
general, which can be independent of the details of anti¬ 
ferrodistortive mode. 

In the point charge model, a semi-quantitative parti¬ 
tion can be applied to the total ferroelectric polarization 
to extract various contributions. The ferroelectric con¬ 
tributions of Ti-O-octahedra, Fe-O-octahedra, and Y-0- 
icosahedra (Fig. |D(a)) can be estimated as^: —144%, 
-1-82%, and -1-162% to the net polarization, as shown in 
Fig. mb). It is clear that the Y-O-icosahedra contribute 
the most polarization, which is partially compensated by 
the Ti-O-octahedra one. The large negative one from 
Ti-O-octahedra can be considered as the induced one 
polarized by the inner-built field from Y-O-icosahedra’s 
ferroelectricity since Ti'^’*' is a proper ferroelectric active 
ion with a considerable large dielectric coefficient. In this 
sense, the ferroelectric polarization in our n = 2 superlat¬ 
tice is not a pure improper one as in the previous studied 
case o^^i^^'^^ , but a hybrid one with both improper and 
proper contributions. 
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FIG. 4. (Color online) (a) Sketch of Y-O-icosahedron (left), 
Fe-O-octahedra (middle), Ti-O-octahedra (right) in the su¬ 
perlattices. (b-c) Point-charge model estimation of individual 
contribution to total polarization at different t/eff’s. Fe^^- 

Ti4+ 

is adopted, (b) Insulating, (c) Metallic. 


D. Possible metallic ferroelectricity 

Very recently, a few studies found some peculiar ma¬ 
terials, in which the ferroelectric distortion persists de¬ 
spite its metallicity^^— . Such a metallic ferroelectricity 
can be characterized according to the structural infor¬ 
mation, e.g. via measuring convergent-beam electron 
diffraction or neutron scattering, as done in Ref. 
Physically, a hint for the metallic ferroelectrics is the 
weak coupling between the electrons at the Fermi level 
and the (soft) phonon(s) responsible for removing inver¬ 
sion symmetry^^— . 

According to the above analysis, the ferroelectric dis¬ 
tortion from Y-O-icosahedra should be robust in our n = 
2 superlattice, which in principle can also be expected 
in other (AR 03 ) 2 /(Ai?’ 03)2 multilayers. Following this 
argument, a promising way to pursuit metallic ferroelec¬ 
tricity is to find a metallic {AB 03 ) 2 /{AB' 0^)2 superlat¬ 
tices. Still taking (YFe 03 ) 2 /(YTi 03)2 as a model sys¬ 
tem, most electron bands near the Fermi level are con¬ 
tributed by Fe^’*' and Ti'^+, but the structural distortions 
of Y-O-icosahedra are the first driving force for polariza¬ 
tion. Thus, it is expectable that the ferroelectric distor¬ 
tion can survive even if the system could become metal¬ 
lic, namely rending a metallic ferroelectricity. 

Despite many experimental approaches (e.g. doping, 
vacancies, or using other perovskites) to make metallic 
{AB 03 ) 2 /{AB' 03 ) 2 , here we use smaller t/eff’s in the 
DFT calculations to obtain the metallicity. Although 
smaller t/eff’s may be not realistic for real YFeOs/YTiOs, 
the propose of this hypothesis is to illustrate the general 
physical mechanism that the improper ferroelectricity in 
(AR 03 ) 2 /(Ai?’ 03)2 is robust against the metallicity, go¬ 
ing beyond a special property limited to a concrete ma¬ 
terial (e.g. stoichiometric YFe 03 /YTi 03 ). In the DFT 


calculation, weak Ues’s give rise to finite density of states 
at the Fermi level although the charge transfer remains 
(almost) complete. Although the Berry phase method 
can not work anymore to evaluate the polarization in the 
metallic state, the point charge model still works. 

With decreasing (C/ti, Upe), the system turns to be 
metallic but its ferroelectric distortion can be even more 
prominent. This enhanced ferroelectricity in metallic 
state is also understandable: the negative contribution 
from Ti-O-octahedra is significantly suppressed while the 
one from Y-O-icosahedra is more robust. Taking the 
(Dxi = 1.6 eV, Upe = 2 eV) case for example, the con¬ 
tribution from Ti-O-octahedra is only 45% of the orig¬ 
inal value, while the Fe-O-octahedra contribution is al¬ 
most unchanged (97% of the original one), as shown in 
Fig. He). The contribution of Y-O-icosahedra is slightly 
increased by 25%, which is the largest contribution to the 
total net polarization. This semi-quantitative partition 
can help to understand the metallic ferroelectricity. 

Our DFT study supports the argument that the im¬ 
proper ferroelectricity due to the geometry factor is ro¬ 
bust against the finite density of states at the Fermi level, 
when the origins of ferroelectric displacement and metal¬ 
licity are different. Although the C/eff’s used here may 
be a little lower for concrete YFe 03 /YTi 03 , the gen¬ 
eral physics raised in the present work is scientific sound. 
Even if the stoichiometric n = 2 YFeOs/YTiOs super¬ 
lattice is not metallic, other approaches can be employed 
to make the superlattice metallic or other systems can be 
designed to realize the metallic ferroelectricity following 
our above argument. 


IV. CONCLUSION 

In summary, the (YFe 03 )„/(YTi 03 )„ superlattices 
have been studied using the standard DFT calculation. 
Since the two parent materials are both Mott insulators, 
unexpected charge transfer has been found, in opposite 
to the intuitional band alignment scenario. In addition, 
the ferroelectricity is predicted in the n = 2 superlat¬ 
tice from the symmetry analysis and confirmed by cal¬ 
culations. Considering the magnetism of Fe, this n = 2 
superlattice is multiferroic. In addition, this ferroelec¬ 
tricity is robust against the metallicity. Even if the real 
(YFe 03 ) 2 /(YTi 03)2 may be insulating, this design prin¬ 
ciple can be extended to other superlattices to search for 
metallic ferroelectrics. 
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